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We calculate the distribution of the scattering matrix at the Fermi level for chaotic normal- 
superconducting systems for the case of arbitrary coupling of the scattering region to the scattering 
channels. The derivation is based on the assumption of uniformly distributed scattering matrices 
at ideal coupling, which holds in the absence of a gap in the quasiparticle excitation spectrum. The 
resulting distribution generalizes the Poisson kernel to the nonstandard symmetry classes introduced 
by Altland and Zirnbauer. We show that unlike the Poisson kernel, our result cannot be obtained by 
combining the maximum entropy principle with the analyticity-ergodicity constraint. As a simple 
application, we calculate the distribution of the conductance for a single-channel chaotic Andreev 
quantum dot in a magnetic field. 
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I. INTRODUCTION 

Statistical aspects of electronic transport through 
chaotic cavities (quantum dots) can be efficiently de- 
scribed using a random matrix model for the N x N 
unitary scattering matrix S of the system (see Ref. [l] for 
a review). For sufficiently low temperatures and voltages, 
transport properties can be expressed by the scattering 
matrix at the Fermi level. Besides unitarity, a crucial 
role is played in the random matrix models by the ad- 
ditional constraints satisfied by S, defining the so-called 
symmetry classes. 

In the absence of superconductivity, following Dyson^ 
three symmetry classes are distinguished, depending on 
the presence or absence of time-reversal and spin-rotation 
symmetry. In this classification scheme, the cases are la- 
beled by the index (3, and the additional constraints on S 
are as follows. In the presence of time-reversal, as well as 
spin-rotation symmetry (f3 = 1), S is symmetric, S — S T . 
In the absence of time-reversal symmetry (f3 = 2), the 
only requirement is the unitarity of S. In the presence 
of time-reversal symmetry, but without spin-rotation in- 
variance (/3 = 4), S is self-dual, S = S R . (The dual of 
a matrix A is defined by A R = tA t t t , with r = i<j2, 
where aj denotes the j-th Pauli matrix in spin space.) It 
has been shown 3 - that for ideal coupling of the scattering 
channels to the cavity, i.e., in the absence of direct re- 
flection from the cavity openings, the distribution of S is 
uniform. Uniformity is understood with respect to the in- 
variant measure in the unitary group subject to the con- 
straints imposed by the symmetries under spin-rotation 
and time-reversal. From the uniform distribution at ideal 
coupling, it follows 3 - that at arbitrary coupling the prob- 
ability density of S is given by the Poisson kernel^^^ 7 ^ 

Pf,(S)<x |det(l-rtS)r^ v+2 -« (1) 

where r is the matrix describing the direct reflections 
from the openings. 

Dyson's classification scheme becomes insufficient 
in the presence of superconductivity^^ In normal- 
superconducting hybrid systems, the scattering ma- 



trix acquires an electron-hole structure ; 11 ! 12 ' 13 and it 
satisfies^ a constraint at the Fermi level, S — Si S* Si, 
expressing the electron-hole symmetry. (Sj denotes the 
j-th Pauli matrix in electron hole space.) Altland and 
Zirnbauer— showed that depending on the symmetries 
under time-reversal and spin-rotation, these systems fall 
into four new symmetry classes, which they labeled fol- 
lowing Cartan's notation of the corresponding symmetric 
spaces. Systems where both symmetries are broken, be- 
long to class D. If only spin-rotation invariance is broken, 
class Z)III is realized. If only time-reversal symmetry is 
broken, the system belongs to class C, and finally, if all 
symmetries are present, the system belongs to class CI. 
The requirements for S following from time-reversal and 
spin-rotation symmetry are the same as in the absence 
of superconductivity. Assuming gapless quasiparticle ex- 
citations, Altland and Zirnbauer introduced a random 
scattering matrix model for transport in chaotic normal- 
superconducting systems, by adopting a uniform distri- 
bution for the scattering matrix. This is appropriate for 
the case when the coupling of the cavity to the transport 
channels is ideal. The analogue of the Poisson kernel, i.e., 
the distribution of S, for the case of arbitrary coupling, 
to the best of our knowledge, has not been presented 
yet. In this study, we aim at providing this result. We 
believe that the knowledge of this distribution is desir- 
able, as it can serve as a starting point to extend results 
that are based on the Poisson kernel PpiS), from Dyson's 
standard symmetry classes to the classes of Altland and 
Zirnbauer. As a particular example we mention the study 
of dephasing in the framework of Biittiker's dephasing 
lead model i 14 ' 15 ' 16 ' 17 In this model, in order to account 
for dephasing mechanisms that occur uniformly in the 
quantum dot, the knowledge of the distribution of the 
scattering matrix for nonideal coupling is essential^ 

The paper is organized as follows. In the next Sec- 
tion, we relate the attributes of the scattering matrix to 
those of normal-superconducting quantum dots, and we 
briefly discuss the conditions at which a random scatter- 
ing matrix description for the transport in such systems 
is adequate. In Sec. IIIIl we detail the properties of the 
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manifolds in the space of N x N matrices defined by 
the constraints on the scattering matrix corresponding 
to the symmetry classes of Altland and Zirnbauer. In 
Sec. IIV1 we present the calculation of the distribution 
P(S) of the scattering matrix, based on the assumption 
that S is uniformly distributed in the ideal coupling case. 
We illustrate the use of our result in Sec. |V]on a simple 
but physically realistic example, a single mode normal- 
superconducting quantum dot in magnetic field. We con- 
clude in Sec. lVII bv contrasting P(S) and the Poisson ker- 
nel Pp(S) regarding the applicability of the analyticity- 
ergodicity constraint of Ref. @. 



II. PHYSICAL REALIZATION OF THE 
SCATTERING MATRIX ENSEMBLES 

In the case of the symmetry classes of Altland and 
Zirnbauer, the role of the chaotic cavity is played by 
a so-called chaotic Andreev quantum dot ft 10 ' 19 i.e., a 
structure formed by a chaotic normal conducting quan- 
tum dot contacted to superconductors. In the vicinity 
of the Fermi level, there are no propagating modes in 
the superconductors. We consider the situation when 
the Andreev quantum dot is contacted to normal reser- 
voirs. The number of propagating modes in the contacts 
to the normal reservoirs (normal contacts for short), in- 
cluding electron-hole degrees of freedom, defines the size 
of the scattering matrix S. We concentrate on the regime 
where transport properties can be expressed in terms of 
the scattering matrix at the Fermi energy: the tempera- 
ture and the voltages applied to the normal reservoirs are 
assumed to be much smaller than the energy scale corre- 
sponding to the escape rate from the normal region and 
the gap of the superconductors. (The superconductors 
are assumed to be grounded.) A sketch of an Andreev 
quantum dot with two superconducting and one normal 
contact is shown in Fig. |TJ (Charge transport can already 
take place using one normal contact, due to the Andreev 
reflection at the superconducting interfaces. 11 ' 12 ' 13 ' 20 ) 

By slightly varying the shape of the Andreev quantum 
dot, one obtains an ensemble of systems and, therefore, 
an ensemble of scattering matrices. We discuss below the 
conditions at which this ensemble can realize the random 
scattering matrix models discussed in this paper. The 
only parameters that enter the scattering matrix distri- 
bution are the symmetries of S and the properties of the 
normal contacts. This implies that the conductance of 
the superconducting contacts should be much larger than 
of the normal contacts, otherwise transport properties 
would be sensitive to the ratio of these conductances! 21 ! 22 
In addition, Frahm et a^ 3 - has shown that for the effect of 
the superconductors on the dynamics in the cavity to be 
considerable, the Andreev conductances of the supercon- 
ducting contacts should be much larger than unity.— For 
a random scattering matrix description of transport, it is 
important that the quasiparticle excitations are gapless. 
If the excitations were gapped, the normal contacts to the 




Figure 1: An Andreev quantum dot formed by a normal con- 
ducting cavity (d) with two superconducting contacts (S). In 
the transport state, an infinitesimal voltage V is applied be- 
tween the grounded superconductors and the normal reservoir 
(N) contacted at the left opening. The distribution of S de- 
rived in this paper considers the effect of a tunnel barrier 
in the normal contact, indicated by a black rectangle in the 
figure. 



Andreev quantum dot would effectively act as normal- 
superconductor interfaces ) 21 ' 22 directly reflecting incom- 
ing quasiparticles, i.e., during transport, the quasiparti- 
cles would not explore the chaotic cavity. Gapless chaotic 
Andreev quantum dots that belong to class C and D can 
be realized with one superconducting contact already, us- 
ing a time reversal breaking magnetic field to suppress 
the proximity gapi 9 ' 23 ' 24 ' 25 For classes CI and Dili, time 
reversal invariance requires the absence of magnetic fields 
in the dot. The gap can be suppresse d 10 ' 25 by using two 
superconducting contacts, with a phase difference n. The 
assumption of a uniformly distributed scattering matrix 
corresponds to assuming that the coupling of the cavity 
to the transport channels is ideal, i.e., that the normal 
contacts are without a tunnel barrier. (The contacts to 
the superconductors can contain tunnel barriers, as long 
as they satisfy the aforementioned requirements for their 
conductances.) In the remaining part of the paper, our 
task is to generalize this uniform distribution to one that 
accounts for nonidcal normal contacts. It is worthwhile 
to note here that we do not rely on the specific details 
of the barriers in the normal contacts, we only use that 
the scattering matrix of the barriers satisfies the same 
symmetry requirements as the scattering matrix of the 
system without the barriers. Our calculation is there- 
fore equally valid for contacts to the normal reservoirs 
with tunnel barriers that do not mix electrons and holes, 
and for barriers that mix electrons and holes. The latter 
situation can occur if there is a region with an induced 
superconducting gap in the contact to a normal reservoir, 
that the quasiparticles have to tunnel through to reach 
the (gapless) cavity region. 
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III. SCATTERING MATRIX MANIFOLDS 

The scattering matrix can be considered as a point of 
a manifold M.x in the space of N x N matrices, where 
X refers to the symmetry class under consideration. The 
distribution of the scattering matrix is understood with 
respect to the invariant measure on A4x- We first state 
the symmetry properties of A4x following Ref. |T(J We 
then take a common route^*^ and consider M.x as a 
Riemannian manifold, to give expressions for the invari- 
ant arclength ds x — Tr(dUdU^) and the correspond- 
ing measure dfj,x(U) in an infinitesimal neighborhood of 
U G .Mx- As usual ) 26 ' ' 28 we parametrize this infinites- 
imal neighborhood with the help of infinitesimal matri- 
ces SUx, with symmetry properties dictated by those of 
Mx, such that the measure is simply the product of the 
independent matrix elements of SUx- 

For class D, the manifold M.d is isomorphic to SO (TV), 
with U= uOu*, O £ SO(iV), and 



1 

v = — 
2 



1+i 1-i 
1-i 1 + i 



(2) 



It might be worthwhile to note here, that solely from 
the unitarity of U and the symmetry U = £iZ7*£i only 
det U = ±1 follows. In Ref. |T(| the manifold A4d was 
identified through the exponentiation of the Bogoliubov- 
de Gennes Hamiltonian, which leads to det U = 1 due 
to the mirror symmetry of the energy levels around zero. 
(The energies are measured relative to the Fermi level.) 
The invariant arclength and measure can be written using 
SU D = O r dO as 



ds 2 D = Tr(6U D 6U%) = 2^{8U D ) 



2 

kl l 



k<J 



dfi D (U) oc ]J(SU D ) M . 

k<l 



(3) 



Note that SU d is antisymmetric due to the orthogonality 
of O. It is seen that ds 2 D [and consequently d^D(U)] is 
invariant under SU D -> WSU D W T , with W € O(JV). 
Such a transformation also preserves the antisymmetry 
of 5U D . 

The manifold A4nm is spanned by U = UU R , 
U € M.D- It is worthwhile to conjugate with 



V 



y/2 V CTl r 



(4) 



and define 6 = V T v^UvV £ SO(iV). The matrix V is 
chosen such that (vV)^ = {vV) R , fr om which it follows 
that U R — vVO R V T v^ . M.DIU is therefore isomorphic 
to the manifold spanned by O = 66 R , 6 € SO(iV). 
Defining SO = (D T dO, the invariant arclength and mea- 



sure can be written in terms of SUdui = SO + SO R as 
ds 2 mn = Tr(SU Dm SUl m ) = Aj2(da 2 kl + db 2 kl ), 



k<l 



d(J,Dia(U) oc W_da kl db kh 



(5) 



k<l 



where, in spin grading 
da db 



SU 



diii 



db —da 



, da = -da T , db = -db T . (6) 



The parametrization ^ follows from SUdui = ~SU]j lu 
and SUdui = SU R U1 . The arclength ds D1II 
and the measure dfiDiu(U) are invariant under 
SUdui -> W5U Dm W T , with W E 0(N). If W also sat- 
isfies W R — W^ 1 , such a transformation preserves the 
symmetries of SUdui- 

In the case of the classes C and CI we omit the spin 
degree of freedom, and we use N to denote the size of 
the scattering matrices without spin. Electron-hole sym- 
metry is now expressed by the relation^ U = E2t/*£2, 
i.e., U is unitary symplectic. For class C this defines 
M.C — Sp(iV). The invariant arclength and measure are 



ds 2 c = Tt(5U c SU f c ) = 



2(EEW ) ) 2 + 2 EEw ) ) 2 )' 



9=1 I 



q=0 k<l 



nn^' nn^ 

\g=l I ) q=0k<l 

Here, SU C = WdU, with 



SU C = l (oh) 5C/ (0) +iJ2^ q su {q \ 

9=1 



(9) 
kl ■ 



(7) 



(8) 



(9) 



where l( ch ) is the identity matrix in electron- hole space, 
and SU^ are N/2 x N/2 dimensional real matrices. Due 
to (SUcY = -SU C : they satisfy <H7(°) = -{SU^) T , 
and for q > 0, 5U {q) = (SU^) T . The arclength and 
the measure are invariant under SUc — * WSUcW^ , with 
W 6 Sp(iV). Such a transformation preserves the sym- 
metries of SUc as well. 

The manifold Aici is spanned by U = UU T , with 
U G Sp(AT). Defining SU — WdU and decomposing it 
according to Eq. ©, we define 

SUci = SU+ {SU) T = zEi(5f7 (1) + iT, 3 SU {3) . (10) 

The invariant arclength and measure are 

ds 2 cl = Tr(SU C iSU CI ) = 



2(EEW) a + 2 EE(*W)« 

9=1,3 I 9=1,3 k<l 



(11) 
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d m (C/)cx J] l[SUif. (12) 

g=l,3 k<l 

The arclength and the measure are invariant under 
SUci -> WflJciW 1 with £ Sp(iV)nO(./V). The sym- 
metry of 5Uci is also preserved under such a transforma- 
tion. 



where we used that in Afs , dfi(So) oc Yljdxj. On 
the other hand, the probability of Afs is the same 
as of Afs , i.e., P(S)d^i(S) = d^i(So), which gives 
P(S) oc | det g(S)\^ 1 ^ 2 , the distribution we are after. 

Parametrizing Afs with the help of SUx and Afs using 
(SUx)o, the relation in Eq. (Til)]) can be written as 



IV. SCATTERING MATRIX DISTRIBUTION 

The scattering matrix S at nonideal coupling can be 
represented as a combination of a random N x N scatter- 
ing matrix So at ideal coupling and a fixed 2N x 2N scat- 
tering matrix S c responsible for the direct reflections. 3 
The matrix So is assumed to be uniformly distributed 
with respect to the invariant measure on A4x- The ma- 
trix S c is given by 



SUx = M'(6Ux) M'\ 



(20) 



S r = 



r a 

t r' 



(13) 



Here the dimension of all submatrices is N x N, and all 
of them carries further structure in electron-hole space, 
and, for classes D and Dili, also in spin space. The 
matrix r describes direct reflection from the contact, r' 
describes reflection back to the cavity from the contact 
and t and t' are the transmission matrices to and from 
the cavity, respectively. The scattering matrices Sq and 
S c have the same symmetries. 

The total scattering matrix S is given by 



S = r + t'S (l~r'S y 1 t, 



and the inverse of the relation is 



S = (t')- 1 (S-r)(l-r^S)- 1 ^. 



(14) 



(15) 



We derive the distribution of S from the uniform distri- 
bution of iSq following a similar logic to the calculations 
in Refs. I29ll30l The starting point of the reasoning is the 
relation 



ss = mss m\ 

where SS = S^dS, SSq = SgdSo and 
M = (1-5V) r 1 . 



(16) 



(17) 



The strategy is to express the arclength in an infinitesi- 
mal neighborhood Afs of S as 



where 



M' = vHlv 

M' = V T viU R M(UfflvV 
M' = U T MUn 



for class D, (21a) 
for class Dili, (21b) 
for class CI, (21c) 



and M' — M for class C. Here the matrices U, Uo are 
used to express S and So for class CI and Dili according 
to Sec. Hill i.e., S = UU V , S = U U$ where y — T and 
U, U € Mc for class CI, and y = R and U, U e Md 
for class DHL The matrix M' satisfies 



M' = M'* 

M' = M'* = tM't t , 

M' = S 2 M'*E 2 

M' = M'* = E 2 M'S 2 



for class D, 
for class Dili, 
for class C, 
for class CI. 



(22a) 
(22b) 
(22c) 
(22d) 



The reality of the matrix elements of M' for class D 
and Dill follows from the fact that the set of matrices 
satisfying A = Sit4*Ei is closed under matrix addition 
multiplication and inversion, and that the combination 
lAAv is real. We show the proof of M 1 = tM't t for 
class DHL Because of vV = t(uV)*t t ' , it is enough to 
show that M = tM*t t , where 



M = U R M(U^ = U^S{1 - 5+r)* _1 S3t/b. (23) 
It is easy to see that 
tM*t t = U R tS*(1 - S T r*)(t*)~ 1 T T S Q H (U^ . (24) 



ds 2 (S) = Tr(dSUS) = Tr(SSUS) = ^g i] {S)dx i dx ]1 Using Eq. and the sclf duality of S and S c we find 



where {dx{\ denotes the set of independent matrix ele- 
ments of SUx in the parametrization of an infinitesimal 
neighborhood Afs of Sq. [Afs is the image of Afs under 
the mapping (|14p .j This way we can relate the measure 
dfj.(So) of Afs to the measure dfx(S) of Afs asi^ 

dfi(S) oc | det g(S)\ 1/2 Y[dx 3 oc | det g(S)\ 1/2 dfi(So), 

(19) 



S? = Tt*(l-S i r*)- 1 T I SM. 



(25) 



Substituting in and using again the self duality of 

S leads to the desired result. The reality of the matrix 
elements of M' for class CI can be proven following anal- 
ogous steps. The relation M' — ^M'*^ for class C and 
CI follows from the closedness of the set of matrices satis- 
fying A = S2j4*S 2 under matrix addition multiplication 
and inversion. 
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Following from properties (]22|) . the matrix M' has a singular value decomposition 

M' = WDW', 

where 



D 


= diaj 


T ,{dk), 


fc = 


1. . 


.iV 


for class D, 


D 


= dia£ 


;(4) i (sp) , 


jfe = 


1 . 


N 
"~2 


for class Dili, 


D 


= dia^ 


;(4) i (ch) , 


k = 


1. 


N 
"~2 


for class C, 


D 


= diaj 


5(4) l (ch) , 


k = 


1. 


N 


for class CI, 



with l( sp ) being the identity matrix in spin space, and 

W,W'eO(N) for class D, 

W = (W R )- X ,W = {W'*)- 1 e 0(N) for class Dili, 

W, W e Sp(iV) for class C, 

W,W e Sp(AT)nO(iV) for class CI. 

Using the decomposition (|26|) . the invariant arclength reads 

ds 2 (S) = Tr(5U x 5U f x ) = Ti {[D(SU' X ) D}[D(SU' X ) D}^} , 

where we used the parametrization (SUx)o — W' (6U x )qW' . From the properties of W in Eq. 
the matrix (6U' x )o has the same symmetries as (SUx)o- It is easily read off that 



(26) 

(27a) 
(27b) 

(27c) 

(27d) 



(28a) 
(28b) 
(28c) 
(28d) 



(29) 

it follows that 



N N 

Vdet g(S) cx J[ d k d t = JJ d^ 1 = \ det M\ N_1 

k<l k 
N/2 N/2 

\/det g(S) oc \[{dkdif = ]J <- 2 = IdetMlf^ 1 

k<l k 

N/2 N/2 N/2 

Vdct g(S) oc JJ (d fc d,) 4 114 = U d l N+2 = |detMr +1 

fc<; j k 

N/2 N/2 

Vdet g(5) oc JJ(dfcdi) 2 = IT < +2 = I det M l f +1 

k<l k 



for class D, 



for class Dili, 



for class C, 



for class CI. 



(30a) 
(30b) 
(30c) 
(30d) 



r 



The distribution of 5 1 is therefore given by 
P(S)oc\det{l~r^S)\- Wt+a) , 



(31) 



where t = 1 in the absence of time reversal invariance 
and t = 2 otherwise, and cr = —1 in the absence of spin 
rotation invariance and a — 1 otherwise. 



V. CONDUCTANCE DISTRIBUTION FOR AN 
ANDREEV QUANTUM DOT IN A MAGNETIC 
FIELD. 

To illustrate the use of our result, we calculate the con- 
ductance distribution for a chaotic Andreev quantum dot 



in a magnetic field. We assume that the spin-orbit scat- 
tering is negligible, i.e., the system belongs to symmetry 
class C. For simplicity, we consider N = 2, which is the 
minimal dimension of S due to the electron-hole struc- 
ture. This corresponds to the case that the quantum dot 
is connected to a normal reservoir via a single mode point 
contact. A sketch of the system is shown in the inset of 
Fig. [21 The point contact is assumed to contain a tunnel 
barrier of transparency T. The barrier alone does not 
mix electrons and holes, therefore its reflection matrix is 
diagonal in electron-hole space, 







) = Vr^Texp(z££ 3 ). 



(32) 



6 



Here £ is the phase an electron acquires upon reflection 
from the barrier. The total scattering matrix S is dis- 
tributed according to P(S) in the group Sp(2) = SU(2). 
The conductance in units of 4e 2 /h is given b y 11 ' 12 ! 13 



G(S) = 1 5, 



he 



(33) 



Writing the total scattering matrix as S = exp(z££3) U, 
U G SU(2), and using that |S he | 2 = \U hc \ 2 and 
d/j,c{S) = d^c{U), the conductance distribution is given 
by 



P(G) 



S(G-\U he \ 2 ) 



-dnc(U). (34) 



/su(2) |det(l-x/T^Tf/)| 3 ' 
Using the Euler angle parameterization for SU(2), 



U= 



e -i{4>+i>)/2 cos (0/ 2 ) - e Ki>-<i>)n s in(6»/2) 
e i(^-V-)/2 sin (g)/2) e l ^+^/ 2 cos(9/2) / 



,i/j,6)e [0, 2tt] x [0, 4tt] x [0, tt] = V, 



(35) 



the measure is d\ic c< sin((9), and |£/he| 2 = sin 2 (0/2). The 
integral 



P(G) 
F r (6,<p,i>) : 



d9d(/)d^Fr(0,(l),ij) 



v 



sin(6l)^[G-sin 2 (6>/2)] 
[2 - T - 2\/r^rcos co7f]3 



(36) 



can be evaluated in closed form, resulting in 

p(n r3 r 2 + 2(G-3)(r-i) 
v ; [r 2 - 4G(r - 1)] 5 / 2 



(37) 



for < G < 1, and otherwise. In Fig. [2] we show P(G) 
for different values of the barrier transparency. It is seen 
that the uniform distribution P(G) = 1 corresponding 
to ideal coupling (r — 1) is gradually transformed into a 
distribution that is peaked at G = as the transparency 
decreases. The first two moments of the conductance are 
given by 



(G) 



r 2 

2 : 



(G 2 



3 ' 



(38) 



VI. CONCLUSION 

In conclusion, we have calculated the distribution P(S) 
of the scattering matrix at the Fermi energy for chaotic 
Andreev quantum dots in the nonstandard symmetry 
classes of Altland and Zirnbauer. Our result, which al- 
lows for arbitrary coupling to the transport channels, is 
based on the assumption that the scattering matrix is 
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G[4e/h] 



Figure 2: Conductance distribution for a single mode chaotic 
Andreev quantum dot in a magnetic field for different values 
of the barrier transparency F. A sketch of the system is shown 
in the inset. The barrier is indicated by a black rectangle. 
The solid, dashed, and dotted curves correspond to Y — 0.9, 
r = 0.6, and T = 0.5, respectively. 



uniformly distributed in Mx for the case of ideal cou- 
pling, i.e., in the absence of direct reflections from the 
openings of the Andreev quantum dot. 

Apart from the symmetry class dependent exponent, 
our result P(S) has a similar structure to the Poisson 
kernel distribution Pp(S) corresponding to Dyson's stan- 
dard symmetry classes. As a closing remark, we would 
like to emphasize an aspect in which P(S) and Pp(S) are 
different. Pp(S) can be obtained^ as the distribution of 
unitary matrices (with the symmetry corresponding to /?) 
that maximizes the information entropy, subject to the 
constraint (S p ) = (S) p = r p , where r is a subunitary ma- 
trix. This analyticity-ergodicity constraint follows 6 from 
the requirement that the scattering matrix has poles only 
in the lower half of the complex-energy plane (analyt- 
icity), and the assumption that spectral averages equal 
ensemble averages (ergodicity). Given the similar form 
of P(S) and P ( a(5), one might wonder whether P(S) can 
be obtained by the same maximization procedure. As we 
show below, the answer is negative. In the presence of 
superconductivity, spectral average can mean two types 
of averages. First, it can refer to averaging scattering ma- 
trices over an interval of excitation energies e. Since elec- 
tron hole-symmetry relates scattering matrices at e and 
—e, it results in the additional constraint S — EiS*!]! 
only at e = 0. Therefore, such an average would be 
over scattering matrices with different symmetry than 
the matrices in the ensemble corresponding to e = 0, 
which violates the ergodicity assumption. Second, the 
spectral average can refer to an average over an inter- 
val of Fermi energies of the superconductor, while the 
excitation energy is kept at s = 0. However, if £ is a 
pole of the scattering matrix as the function of the Fermi 
energy, so is £*. To see this, one turns to the channel 
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coupled model used in Ref. 

IE 

in which the poles of the 
scattering matrix on the complex-Fermi energy plane are 
eigenvalues of a matrix (ft — iWW^Ss, where ft mod- 
els the Bogoliubov de Gennes Hamiltonian (at a fixed 
Fermi energy) and W is a coupling matrix. The conju- 
gation relation between the poles is the consequence of 
Ei (ft - iWW^)*T,i = -(ft - iWW+). This precludes 
the use of the analyticity properties as in Ref. 0. We 
thus find that the analyticity-ergodicity constraint is not 
applicable to the ensembles studied in this paper. This 
is already signaled in the case of ideal coupling. For ex- 
ample for class D, one has±£ (S) = 0, but (TrS 2 ) = 1. 
For the more general, nonideal coupling case described 



by P(S), one also finds that r ^ (S), in contrast to 
the case of the Poisson kernel Pp(S). This can be il- 
lustrated on the example in Sec. El where, for £ = 0, 
(TrS) = VT" 7 T(r + 2), as opposed to Tr r = 2\/l — T. 
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